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Abstract 

The thermal evolution of a neutron star is studied by including the energy loss due to axion emission. Two 
axion models and three types of neutron-star matter equation of state are used with the effects of nucleon 
superfluidity properly taken into account. In comparison with the observational data of PSR0656+14 from 
ROSAT, the upper limits on the axion mass are found to be m a < 0.06 — 0.3 eV and 0.08— 0.8 eV for 
the KSVZ and DFSZ axion models, respectively, with the soft equation of state giving the most stringent 
limits. 



The axion arises as a solution to the strong CP problem 
(Turner 1990, Raffelt 1990). While the standard axion 
model was excluded by experiments, the invisible axion 
model has survived mainly because the axion's coupling 
to matter is weak, which is an unknown parameter in the 
theory. Over the years, various laboratory experiments 
as well as astrophysical arguments have been used to con- 
strain its parameters. Since laboratory experiments can 
explore only a limited parameter regime, including those 
planned in the foreseeable future, astrophysical consider- 
ations have played an important role in placing the limits 
on the axion parameters. Within these limits, the axion 
remains as one of the candidates for dark matter. 

There arc two types of axion models — the KSVZ 
(hadronic) model (Kim 1979, Shifman et al. 1980) and 
the DFSZ model (Dine et al. 1981, Zhitnitskii 1980). In 
the KSVZ model, the axion couples only to the photons 
and hadrons, while in the DFSZ model the axion cou- 
ples to the charged leptons as well. The axion-fermion 
and axion-photon coupling constants as well as the axion 
mass are unknown parameters in these theories. Cur- 
rently, cosmological arguments give m a > 10~ 5 eV (Ab- 
bott and Sikivic 1983, Dine and Fischler 1983). The limit 
from Supernova 1987A, which used to give m a < 10~ 3 
eV, is now somewhat relaxed m a < 0.01 eV (Raffelt and 
Scckel 1991, Janka et al. 1996). The red giant limit 
m a < 0.009/ cos 2 (3 (Raffelt and Weiss 1995) applies only 
to the DFSZ model. The laboratory experiments give 



weaker limits. 

In the present paper we study how axion emission af- 
fects the thermal evolution of neutron stars. We use 
the neutron star evolutionary code with three types of 
equation of state to calculate the surface temperature 
of neutron stars. We compare theoretical cooling curves 
with observation and obtain the upper limits on the ax- 
ion mass, which are weaker than, but comparable with, 
the limit from SN 1987A. 

Axion Emissivity: In neutron stars, the dom- 
inant axion emission mechanisms are the following 
bremsstrahlung processes in the stellar core: n + n — > 
n + n + a, p + p — > p + p + a, and n + p^n + p + a, 
where n,p, and a are the neutron, proton and axion. The 
energy loss rate of each process, in the units Ti = c = 1, 
is given by (Iwamoto et al. 1998), 
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9 = 9ap P + gann, h = g app - g ann ; x = m 7T /2p F (n), y = 
/2p F (p); f ~ 1 is the pion-nucleon coupling constant; 
p F (n) ~ 340(p/p ) 1/3 MeV/c, p F (p) ~ 85(p/p ) 2/3 
MeV/c are the nucleon Fermi momenta; m p , m n and m T 
arc the proton and neutron effective masses and pion 
mass, respectively. 
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is the axion-nucleon coupling constant, where i = p (pro- 
ton) or n (neutron), m N is the nucleon mass, and f a is 
the axion decay constant. Cj depends on the models: the 
DFSZ model gives 

c p = -0.10 -0.45 cos 2 (i, c n = -0.18 + 0.39 cos 2 j3, (7) 

and the KSVZ (hadronic) model gives 

Cp = -0.385, c„ = -0.044. (8) 

The axion mass is related to f a via 
0.0074 
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We note that axion emission is suppressed if nucleons 
become superfluid, as in the case of neutrino emission 
involving nucleons. 

Results and Discussion: We employ the numerical 
calculation code essentially the same as the one described 
in Umeda et al. (1994), except for the inclusion of the 
energy loss due to axion emission. We neglect internal 
and other possible heating mechanisms as well as the ex- 
istence of non-standard cooling mechanisms. The baryon 
mass of the neutron star is set to 1.4 M©. 

Theoretical cooling curves are compared with the 
observational data for three pulsars: PSR 1055-52, 
Geminga and PSR 0656+14 (sec Tsuruta 1998 and 
Becker 1994 for references). The energy loss rate due to 
axion emission is proportional to the axion mass squared, 
rv? a \ therefore, we can obtain the upper limit on the axion 
mass from the condition that the cooling curve does not 
pass below the lower bounds on the observational points. 
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Fig. 1. Cooling curves with KSVZ axion emission 



In Figure 1, we show the standard cooling curve and 
those with the KSVZ axion model for four different ax- 
ion masses (or f a ). The FP equation of state and the TT 
neutron 3 P2 superfluid energy gap (Takatsuka and Tam- 
agaki 1993) are adopted. Since the data point for PSR 
1055-52 is located above the standard cooling curve, we 
do not use this data: this is likely to be due to some other 
(unknown) effects. Conservative limits can be obtained 
by using the other two data. Figure 1 shows that the 
PSR 0656+14 gives a more stringent limit than Geminga, 
and hence we obtain the axion mass limit from the lower 
bound on the PSR 0656+14 data. 

The results for both the KSVZ and DFSZ axion models 
with stiff (PS), medium (FP) and soft (BPS) equations 
of state are summarized in Figures 2-4. The BPS model 
gives the more stringent limit. This is because the TT 
gap vanishes in the high density region (i.e., inside the 
stellar core) with this equation of state; thus, axion emis- 
sion is not suppressed. Extending superfluidity to higher 
density regions will have an effect similar to increasing 
the stiffness of the equation of state. For example, in 
the FP model, if the AO neutron 3 P 2 gap (Amundsen 
and 0stgaard 1985) is adopted, m™ ax is 0.3 eV, while if 
there is no neutron 3 P 2 superfluid, m™ ax is 0.06 eV. Note, 
however, that the AO model probably overestimates the 
energy gap at high densities, because the density depen- 
dence of the neutron effective mass is neglected. Future 
refinements of the observation will provide more stringent 
limits. 
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Fig. 2. Axion mass limits (m™ ax ) in the FP model 




Fig. 3. Same as Fig. 2 in the PS model 
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Fig. 4. Same as Fig. 2 in the BPS model 



